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We present the first experimental tests of the recently derived fusion rules for Navier-Stokes (N-S) 
turbulence. The fusion rules address the asymptotic properties of many-point correlation functions as 
some of the coordinates coalesce, and form an important ingredient of the nonperturbative statistical 
theory of turbulence. Here we test the fusion rules when the spatial separations lie within the inertial 
range, and find good agreement between experiment and theory. An unexpected result is a simple 
linear law for the Laplacian of the velocity fluctuation conditioned on velocity increments across 
large separations. 


In this Letter the predictions of the recently proposed 
fusion rules [Q| are tested by analyzing the turbulent ve¬ 
locity signal at high Reynolds number. We start with a 
short theoretical summary. 

The theory focuses on two-point velocity differences 

w(r , r', t) = u(r', t) — u(r , f), (1) 

where u(r' ,t) is the Eulerian velocity field accessible to 
experiment. One attempts to extract predictable and 
computable results by considering the statistical proper¬ 
ties of w H|jj . The most informative statistical quantities 
are the equal-time rank-n tensor correlation functions of 
velocity differences 

Fn(r 1 ,r' 1 -,r 2 ,r' 2 -,.. .;r n ,r’ n ) 

= (w(r 1 ,r' 1 )w(r 2 , r' 2 )... w(r n , r’ n )) , (2) 

where (•) denotes averaging, and all coordinates are dis¬ 
tinct. We consider N-S turbulence for which the scaling 
exponents are presumed to be universal (i.e., they do not 
depend on the detailed form of forcing), and where the 
correlations IF n are homogeneous functions Q , namely 

^n(Ari, Xr [,..., \r' n ) = \ Cn J r n (r 1 ,r ' 1 ,..., r' n ), (3) 

( n being the homogeneity (or scaling) exponent. This 
form applies whenever all the distances | rq — r'| are in 
the “inertial range”, i.e. between the outer scale L and 
the dissipative scale ry of the system. Our aim here is 
to describe the behavior of such functions as pairs of co¬ 
ordinates approach one another, or “fuse”. The fusion 
rules, derived in |]I],Q, govern the analytical structure of 
the correlation functions under this coalescence. 

The statistical function that has been most commonly 
studied [^]|3||||]] is the structure function S n (R) 

S n {R) = (| w(r, r ')\ n ), R = r'-r. (4) 

Clearly the structure function is obtained from (@) by 
the fusing of all coordinates rq into one point r, and all 


coordinates r\ into r' = r + R. In doing so, one crosses 
the viscous dissipation length-scale. One then expects a 
change of behavior, reflecting the role of the viscosity in 
the theory for S n (R). In developing a N-S based theory 
in terms of S n (R), one encounters the notorious closure 
problem: one must balance terms arising from the con¬ 
vective term it ■ Vit and the dissipative vYJ 2 u term, nei¬ 
ther of which can be neglected. Hence determining S n (R) 
requires information about S n +i(R). All known closures 
of this hierarchy of equations are arbitrary. However, ac¬ 
cording to the theory of Refs. QQ, the fully unfused T n 
does not suffer from this problem. When all separations 
are in the inertial range, the viscous term may be ne¬ 
glected, and one obtains [Q homogeneous equations for 
T n in terms of fF n only, with no hierarchic connections 
to higher or lower order correlations. Such homogeneous 
equations may exhibit new, anomalous scaling solutions 
for the correlation functions fF n . 

There are various possible configurations of coales¬ 
cence. We will test only those fusions in which the coa¬ 
lescing points are those of velocity differences. One can 
also consider the coalescence of points from different ve¬ 
locity differences, but they are experimentally more dif¬ 
ficult to measure; we will “precoalesce” all such points 
here and comment on its effects Q]. The first set of fu¬ 
sion rules that we examine concerns IF n when p pairs 
of coordinates n, r[ ... r p , r' p , ( p < n) of p velocity dif¬ 
ferences coalesce, with typical separations between the 
coordinates r,; — r' ~ r for i < p, and all other sepa¬ 
rations of the order of R, r <C R <C L. The fusion rules 
predict 

R n (r 1 ,r' 1 ;...;r n ,r , n ) (5) 

= Fpiri, r i;...; r p , r' p )^! n . p {r p+l , r' p+1 - ...;r n , r'J , 

where IF p is a tensor of rank p associated with the first p 
tensor indices of IF n , and it has a homogeneity exponent 
C p- The (n - p)-rank tensor '& ntP (r p+1 ,r p+1 ;... ; r n , r'J 
is a homogeneous function with a scaling exponent 
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FIG. 1. Log-log plot of the structure functions S n (R) as 
a function of R for n = 2,4,6, 8 denoted by +, x, * and o 
respectively. 

Cn~ Cpi and is associated with the remaining n— p indices 
of J- n . The case p = 1 is special, since the leading order 
evaluation cancels by symmetry. The next-order result, 
for a randomly oriented set of pairs with the separation 
R, is 

Fn(ri, r'i,... ; r n , r' n ) ~ (r/R)S n (R). (6) 

One can also consider correlation functions under the 
operation of gradients. One such set of functions, de¬ 
noted J n (R), have particular significance. They arise 
from the dissipative term when one obtains a statistical 
balance from the N-S equations. For comparison with 
experiments we define J„(i?) in terms of the longitudinal 
velocity difference Sur = w(r, r + R) ■ R/R as 

J n {R) = v(v 2 u{r)[5u R ] n - 1 ). (7) 

The Laplacian operator in (|7j) is interpreted as a finite 
difference of longitudinal components of the velocity, 

V 2 tt(r) = [w{r, r + p) — w(r, r — p)] ■ p/p 3 . (8) 

The predictions of the fusion rules are different for p 
above and below the dissipative scale. They read: 

Jn(R) = nC n J 2 S n {R)/2S 2 {R) , p > V , (9) 

J n {R) = nC n J 2 S n+1 (R)/S 3 (R) , p <C 77 , (10) 

where C n and C n are R-independent dimensionless con¬ 
stants. The fusion rules do not rule out an n-dependence 
of these coefficients. J 2 is expected to be R-independent. 

Here we will test these predictions using atmospheric 
turbulence data obtained by means of a single hot-wire 
probe mounted at a height of 35 m on the meteorological 
tower at the Brookhaven National Laboratory. The hot¬ 
wire was about 0.7 mm in length and 6 pm in diameter. 
It was operated on a DISA 55M01 anemometer in the 
constant temperature mode. The wind direction, mea¬ 
sured independently by a vane anemometer, was approx¬ 
imately constant. The frequency response of the hot-wire 
was good up to 5kHz. The voltage from the anemometer 
were low-pass filtered at 2 kHz and sampled at 5 kHz, and 
later converted to velocity through an in-situ calibration. 
The mean wind speed was 7.6 ms -1 and the root-mean- 
sqaure velocity was 1.3 ms -1 . With the usual procedure 
of surrogating time for space (“Taylor’s hypothesis”), we 
obtain the Taylor microscale Reynolds number to be 9540 
and the Kolmogorov microscale to be 0.57 mm. 

In Fig. [l]we present the structure functions S n (R ) as a 
function of R. They were computed using 10 million data 
samples. In this and other figures, spatial separations 
have units of sampling times, and the velocity is nor¬ 
malised by the RMS velocity. This figure shows that we 


have around three decades of “inertial range” (between, 
say 10 and 10 4 sampling units) but that the dissipative 
range is not well-resolved. Structure functions of orders 
higher than 8 are less reliable and will not be considered. 

While the fusion rules are formulated for differences in 
d-dimensional space, the surrogated data represent a 1- 
dimensional cut. This has implications for the choice of 
positioning of the coordinates. In d-dimensional space we 
can choose separations to fall within balls of size R and 
r respectively. In our case this ball collapses onto a line, 
and best results are obtained when the pairs of coordi¬ 
nates in the two groups coincide. As a simple demonstra¬ 
tion consider the second order quantity J r 2 {r\,r' 1 ,r 2 , r 2 ) 
with the three different choices: (i) r\ = r[ = x, 
r 2 = r ' 2 = y, r = \y - x\, (ii) |ri - r[\ = \r 2 - r' 2 \ = r, 
r[ = r 2 , (iii) |ri — r(| = \r 2 — r' 2 \ = r where r[ and r 2 
are also separated by r. In one dimension one can simply 
compute the ratio of the correlation functions in cases (ii) 
and (iii) with respect to case (i). One finds a reduction 
factor 2^~ l -l « -0.2 and (-2« 2+1 + l+3 C2 )/2 « -0.05. 
We thus see that there is a rapid decrease in amplitude 
when the distances are not enmeshed, and so all aver¬ 
aging is done using maximally enmeshed configurations 
(i.e., case (i)). As remarked above, in doing so, all the 
“unprimed” points not across velocity differences are al¬ 
ready fused. This procedure does not affect predictions 
(1) and (§), see (§. 

Explicitly, therefore, we examine the behavior of the 
correlation function 

Tp+qfr, R) = ([u(x + r)-u(x)] p [u(x + R) — u(x)] q ) (11) 

as a function of both r and R for several values of the 
powers p and q. From Eq. (||) one expects 

•Wg R) ~ S p (r)S q+p (R)/S p (R). (12) 

FIG. 2. Log-log plot of J-p+ q (r,R) as a function of r at 
fixed R for q = 2 and p = 2, 4, 6 denoted by +, x and * re¬ 
spectively with dashed lines. Shown with dotted lines are the 
same quantities divided by S p (r). 

FIG. 3. As in Fig. 2 as a function of R at fixed r, 
with the dotted lines representing the quantities divided by 
S p+q (R)/S P (R). 

FIG. 4. Log-log plot of T\+ q (r,R) as a function of R for 
q = 1,3,5 denoted by +, x and * respectively. Shown with 
dotted lines are the same quantities divided by Si+ q (R)/R. 


FIG. 5. log 10 [J n (i?)] as a function of the fusion rule pre¬ 
diction \og 10 [J 2 Sn{R)/S 2 (R)] for n = 2,4, 6, 8 denoted by +, 
x, *, o respectively. Inset: the coefficient C n with the same 
notation. 
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In Fig. |]we display the results for q = 2 with even values 
of p as a function of r for r in the inertial range. Only 
even values are displayed as the odd correlations fluctu¬ 
ate in sign. The large scale R was fixed at the upper 
end of the inertial range. The data show clean scaling 
in the inertial range. Overlaid are the averages corrected 
by the prediction of the fusion rule (|l^) . Here and in all 
other figures the averages themeslves are connected with 
dashed lines, whereas compensated results are shown dot¬ 
ted. One observes a change of behavior as r approaches 
R at the upper limits and the average increases in size 
towards the “fully-fused” quantity, S p + q (R). Similarly 
convincing results were obtained for other values of p 
and q. 

In Fig. H we show fF p+q (r, R) as a function of the large 
scale R with the small separation r fixed at the low end 
of the inertial range, together with the values corrected 
by ©■ There is a clear trend towards zero slope in the 
corrected quantity in the upper inertial range. 

We consider now the special case that a single pair 
of points in a velocity difference approach one another. 
The prediction of Eq. d) is tested for fixed r and the 
expected dependence on R is well-verified in Fig. [|. The 
results found by varying r are not shown here: the linear 
configuration of all the measurement points leads to a 
competition between the leading and next order of scal¬ 
ing. 

The function J 2 was computed and confirmed to be 
constant throughout the inertial range. In Fig. |] we 
present J n (R) as a function of nJ 2 S n {R)/2S 2 {R) for n, 
n = 2,4, 6,8 and inertial range R. The finite difference 
surrogate of the Laplacian (cf. Eq. (|j)) was computed 
with p = 10. The straight line y = x passing through the 
data is not a fit. It appears from these results that <D is 
obeyed well with C n = 1. The R independence of C n is 

FIG. 6. The conditional averages of the Laplacian 

a direct confirmation of the fusion rules for the fusion of 
two points. On the other hand, the fact that C n is inde¬ 
pendent of n is a surprise that does not follow from fusion 
rules, and has interesting implications for the statistical 
theory of turbulence. A more sensitive check of the value 
of C n is obtained by dividing nJ 2 S n {R)/S 2 {R) by J n (R) 
for individual values of n and R. This is displayed in the 
inset in Fig. 5. Clearly, there are statistical fluctuations 
that increase with increasing n but the data show that 
C n is approximately constant in R and n, with a value 
of about unity. 

If indeed C n is n- independent, there are surprising con¬ 
sequences for the conditional statistics of our field. To see 
this, J„ (R) may be rewritten 

J n {R) = J d5uRP[SuR}(\/' 2 3 4 5 u(r)\SuR)Su 1 jf 1 . (13) 

Here ( V 2 u(r)\8v,R ) is the average of the finite difference 


Laplacian conditioned on a value of Sur. The only way 
of satisfying both (|ui]) and (§ with C n that is indepen¬ 
dent of n and R is to assert that the conditional average, 
which is in general a function of Sur and R , can be fac¬ 
tored into a function of R and a linear function of 5ur: 

( V 2 u(r)\6u R ) = 2|g ^ Sur. (14) 


Such linear laws have been discussed in the context of 
conditional statistics of passive scalar advection jf] [TT| , 
and were thought to be reasonable because of the linear 
nature of the advection-diffusion equation for the scalar. 
Thus, linear laws for N-S turbulence as well were un¬ 
expected. In Fig. 6 we display a direct calculation of 
the conditional average of the surrogate Laplacian with 
p = 10, multiplied by 2S 2 (R)/J 2 as a function of Su for 
nine values of R spanning the inertial range. All the data 
collapse on a single straight line whose slope is unity. For 
ease of presentation we have displaced the different val¬ 
ues of R from each other. It appears that the prediction 
& is amply supported by the data. 

Unfortunately we cannot test © with the present 
data because sub-dissipation scales are not resolved. The 
prediction implies that the nature of the conditional av¬ 
erage changes qualitatively when p decreases below the 
dissipative scale. Such changes have important conse¬ 
quences for the ultraviolet properties of the statistical 
theory of turbulence, and a rich variety of predictions are 
already available 1121. It is thus worthwhile generating 
high-Reynolds-number data that resolve sub-dissipative 
scales. Efforts to acquire such data are under way. 
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